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Abstract 



Some time ago, an accurate phenomenological approach, the BSW model, 
was developed for proton-proton and antiproton-proton elastic scattering 
cross sections at center-of-mass energies above 10 GeV. This model has been 
used to give successful theoretical predictions for these processes, at succes- 
sive collider energies. 

The BSW model involves a combination of integrals that, while com- 
putable numerically at fairly high energies, require some mathematical anal- 
ysis to reveal the high-energy asymptotic behavior. In this paper we present 
a high-energy asymptotic representation of the scattering amplitude at mod- 
erate momentum transfer, for the leading order in an expansion parameter 
closely related to the logarithm of the center-of-mass energy. 

The fact that the expansion parameter goes as the logarithm of the energy 
means that the asymptotic behavior is accurate only for energies greatly 
beyond any foreseeable experiment. However, we compare the asymptotic 
representation against the numerically calculated model for energies in a less 
extreme region of energy. The asymptotic representation is given by a simple 
formula which, in particular, exhibits the oscillations of the differential cross 
section with momentum transfer. We also compare the BSW asymptotic 
behavior with the Singh-Roy unitarity upper bound for the diffraction peak. 
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1 Introduction 



Forty years ago, it was found theoretically, on the basis of quantum field 
theory, that, contrary to the general belief at the time, the total cross sections 
for hadronic scattering increases monotonically without limit at high energies 
[1] . In order to make predictions that could be verified by later experiments, 
it was essential to develop an accurate phenomenology with the following 
characteristics: 

(i) - it agrees with the above theoretical asymptotic result at high energies, 

(ii) - it describes the experimental data at energies available at that time. 
Such a phenomenological model, termed the BSW model, was formulated 

by three of us El S] . 

In this paper, after reviewing the BSW phenomenological model for the 
elastic scattering amplitude, we present an asymptotic representation for the 
scattering amplitude for extremely high center-of-mass energies. 

In section 2 we recall the basic features of the impact picture approach 
(BSW), with its phenomenological parameters, and the scattering amplitude 
is rewritten in a form appropriate to get its asymptotic expansion. Another 
expression is obtained in section 3, by means of the exact evaluation of the 
function F that occurs in the model as a Fourier transform integral. In 
section 4 we obtain the high-energy asymptotic representation of the BSW 
amplitude in terms of a simple expression. This derivation involves a num- 
ber of mathematical steps and the consideration of two different kinematic 
regions. Some technical details are collected in the Appendix. Numerical 
calculations are presented in section 5, where we discussed some features of 
the real and imaginary parts of the scattering amplitude, generating some 
oscillations in the differential cross section. The comparison of the asymp- 
totic representation and the exact BSW result is done at the LHC energy 
and also at a much higher center-of-mass energy of 6000TeV. In section 6, 
we also compare the asymptotic representation with the Singh- Roy unitarity 
upper bound for the diffration peak and we give our concluding remarks. 
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2 The BSW model 

To describe the experimental data taken at the relatively low energies avail- 
able to experiments forty years ago, the BSW model was proposed, including 
Regge backgrounds. Both for the energies of the present-day colliders and 
for the purpose of studying the asymptotic behavior of the model at high 
energies, all the Regge backgrounds can be neglected. The BSW model is 
given by the following matrix element for elastic scattering 



M{s,A) = —j dx^e'^^-'^^Dis, x^) , (1) 

where s is the square of the center-of-mass energy, A is the momentum 
transfer, x±^ is the impact parameter and all spin variables have been omitted. 
For this model we take the simplest form that we can use for the opacity 

D(s,a;x) = 1 -e-^(^'^^) , (2) 

with 

n{s,x^)=S{s)F{x^) , (3) 

where x± = \x±\ . The function S{s) is given by the complex symmetric 
expression, obtained from the high energy behavior of quantum field theory 

m 

c c 

S{s) = + -\j , (4) 

with s and u in units of GeV^, where u is the third Mandelstam variable. In 
this Eq. (|4]), c and d are two dimensionless constants given below in Table 
1. That they are constants implies that the Pomeron is a fixed Regge cut 
rather than a Regge pole. For the asymptotic behavior at high energy and 
modest momentum transfers, we have to a good approximation 

In M = In s — ivr , (5) 

so that 

17VC 

^^'^ = (M^ + (In.-ZTr)^' ■ 

Because F depends on x± only through x±, the Fourier transform in Eq. ([T]) 
simplifies to 

/•oo 

M{s, A) = is / dx^ XX Mx^A) [1 - e-'SW^(^i)] , (7) 
Jo 
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c 


= 0.167, 


c' 


= 0.748 


nil 


= 0.577 GeV, 




= 1.719 GeV 


a 


= 1.858 GeV, 


/ 


= 6.971 GeV-2 



Table 1: Parameters of the BSW model 



where A = |A|. The function F{x±) is taken to be related to the electro- 
magnetic form factor G{t) of the proton, where t = — is the Mandelstam 
variable for the square of the momentum transfer. Specifically, F{x±) is 
defined as in [2] via its Fourier transform F{t) by 

m = fm)?^^ , (8) 

with ^ 

^^^^ = (l-t/mf)(l-t/mi) • 

The remaining four parameters of the model, /, a, mi and m2, are given in 
Table 1. 

The task is to study the asymptotic behavior of Ai for large In s and modest 
momentum transfers. 



3 The evaluation of F{x±) and its consequences 

The purpose of this section is to find the exact expression of F{x±), entering 
in Eq. ([7]), in order to determine the most relevant region in x± for the 
calculation of the asymptotic limit of A4{s, A), for large s. 
Noting that F depends only on A^, the Fourier transform that defines F 
simplifies to an integral over one variable, so that we have 

POO 

= / dAAF{-A^)Mx^A) , (10) 
Jo 

where Jo denotes the Bessel function of zero order. From Eq. ([8]), we have 
explicitely 

^ ^ ^ (1 + AVm2)2(l + AVm2)2a2 + A2 ' ^ ' 
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which is a rational fraction, symmetric in mi,m2, whose decomposition into 
simple elements, allows the direct calculation of F{x±). As expected, the 
final result can be expressed in terms of modified Bessel functions Kq and 
Ki. We have the decomposition 



(12) 



where 



Fi(x±) 



fm\m\ 



2{m\ — m\Y \ m 



mf + m? 



0? — m\ 



\{miXi_)Ki{miXi_) 



-4 



+ 



2\2 



Ko{miX±) 



(13) 



F2{X^) 



£ 4: 4 

J 77117712 



2(m2 



771 



2\3 



77ll 



mf a + 7712 



+4 



+ 77il 
a? — rrio 



7717, 



[77lo 



7717, 



m2X_L)Ki{77l2X±) 



Ko{77l2X±) 



4^2 



F3{X^) 



2f 77117712a 



mf)^(a^ 



Ko{axA_) 



(14) 
(15) 



which is clearly symmetric in mi, m2- The arguments of the modified Bessel 
functions Ki and Kq are mix^ for Fi, while m2X^ for F2 and ax_L for F3. The 
function F is real, positive, bounded, and monotonically decreasing toward 
as x^ increases without bound. Because RciS is large and positive for large 
values of Ins, the only contribution to Ai comes from x± fairly near Ins, 
so that the large argument asymptotic expressions for the modified Bessel 
functions are applicable. Since as seen from Table 1, mi < m2 < a, F2 and 
-F3 are exponentially smaller than Fi, only Fi contributes to the asymptotic 
behavior of A^(s,A). For the moment we keep the exact expression for 
F(x^), but we change variables to emphasize the role of Fi as follows. Let 
us define two dimensionless variables 



X 

a 



7nix±, 
A/mi. 



(16) 
(17) 



With this change of variables, Eq. ([7]) becomes 



JO 



dxx Jo(ax) 



-S{s)F{x) 



(18) 
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where F{x) = -Fi(x) + F2{x) + Fsi^x) and Fj{x) = Fj{x/mi) = Fj{x_i), 
for j = 1,2,3. We note that for x ^ 1 the asymptotic representation of 
F, imphed by Eqs. f|T3l) and f|T6|) and the asymptotic representation of the 
modified Bessel function [5j, is 



2{ml — mf)2(a2 — m() 
fev^e-^ (20) 



^(^) - ^^i2'lL^^,_., ^K^ix) (19) 



where b is the real coefficient 



fm\m\{a^ + mf) 
2(m2 — m\Y{a? — mf) 



(21) 



4 High-energy asymptotic behavior 

Some additional approximations make it possible to obtain an expression 
of the scattering amplitude in the high energy limit. The integrand of the 
integral in Eq. ([7]) is complex. From the definition of S{s) in Eq. (jl]) and 
the explicit expression in Eq. ([6]), we notice that for sufficiently high values 
of s, namely 

In s > TT , (22) 

the second denominator in Eq. ([6]) can be approximated by the first denom- 
inator, so that 

^(''-(llF-(f)-''"- 

Thus the phase of S{s) is found to be a constant, namely, — |vrc ~ —0.26 
and when In s is large one has 

Re5(s) > Im5 > 1. (24) 

The function F{x) is positive for all real x and decreasing exponentially with 
large x as stated in Eq. (12Ui) . Because of these properties of F{x) and the 
large value of Re5, the opacity is essentially 1 for values of x, well below a 
transition value, and drops exponentially to zero as x increases well above this 
transition value. From this, it follows that the asymptotic representation of 
the scattering amplitude can be obtained by replacing F{x) by F(zo)e~*^^'~^"\ 
where Zq is a complex transition value chosen in such a way to make S{s)F{zq) 
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real and of order 1. 

By means of a translation in the complex plane from x to x' = x — zq, the 
scattering amplitude of Eq. ffTSjl becomes asymptotically 



IS 



m 



M{s,A) ~ — 2 / dx {x' + zq) Jo[a{x' + zo)] 



IS 



1 J -ZQ 





1-e 



-S{s)F{zo)e'- 



dx'{x' + zo)Jo[a{x' + zq)] + A{s, zq) 



1 



2 



where we define 



— Ji(azo) + ^(s,2o) 
a 



dx'{x' + Zo)<>^o[a(a;' + zq)] 



(25) 
(26) 
(27) 



1 _ ^-S{s)F{zo)e- 



dx'{x' + zo)Jo[a{x' + ZQ)]e 



-5(s)F(2o)e- 



(28) 



It remains to determine zo and then to evaluate the asymptotic representa- 
tion of A[s,zo). In order to find zq we use the following relation derived in 
Appendix A (see Eq. flA.7P ) 







dx 



1-e 



dx e 



which suggests defining zq as solution of the equation 

S{s)F{z^) = e-^ , 



(29) 



(30) 



where 7 ~ 0.5772 is the Euler's constant. While zq, as a function of s, is best 
obtained by solving Eq. f l30|) numerically, one can see its approximate value 
by using the asymptotic representation for F, so that 



By taking logarithms one finds 

Zq ^ \n[h^S{s)] +7, 



(31) 



(32) 



showing how Re 2:0 grows with Re ln[5(s)] and how Im^o = arg(5(s)) ap- 
proaches —0.26, as Ins becomes very large (see Fig. 1). Then the fact that 
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Re^o is large can be used twice. First, in Eq. (125]) . one can safely replace the 
lower limit —Zq in one of the integrals by — oo, so that 



A{s,Zo) 



dx{x + Zo)Jo[a{x + Zq)] 



l-e 



dx{x + zo)Jo[a{x + zo)]e~ 



(33) 



Secondly, since the important regions of integration are where x is of the order 
of 1, to determine the leading order behavior of the scattering amplitude it 
suffices to replace the factors {x + zq) by Zq, so Eq. (155]) reduces to 



A{s,zo) 



oo 

~ Zo-{ I dx Jo[a{x + Zq)] 







dx jQ[a{x + Zo)]e 



1-e 



_g-{a: + 7) 



(34) 



To determine the asymptotic behavior of A{s,zq), we have to study the 
Bessel function which depends upon the parameter a, and there are two re- 
gions to consider. 



i) Small a region 

By assuming that a is small, but not azo, we keep only the first three terms 
of the Taylor series to obtain: 

(qx)'^ 

Jo[a{zo + x)] = Jo{azo) — axJi{azo) 



-Jo{azo) + higher order terms. 

(35) 

By substituting Eq. f l35|) into Eq. f p4|) . one obtains after some integration by 
parts 



A{s,zo) 



a. 



2^i(e '^)Jo{azo) - -l2{e '^)Ji{azo) 



a 



—Isie '^)Jo{azo) 





a 



^2 Jii.(^Zo) + yC(3) Jo(tt2;o) 



(36) 



where the X„ are defined in Appendix A and Q is the Riemann Zeta function, 
and therefore C(3) ^ 1.2021. 
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ii) Large a region 

In this region the large argument asymptotic expansion of the Bessel function 
[5] allows one to write, for large \zo\ and the contributing values of x which 
are 0(1), 



Jo[a{zo + x)] = \l cos[a{zQ + x — vr/4)] + higher order terms 

V TT^o 

\/1txzq 

By substituting Eq. fl37|) into Eq. flM|) . this yields for this region of a 




A{s,z^) ~ re*(°"»-"/^)j^(2«) + e-*(""«-"/^)j^(-m)l (38) 

>' 27r 



zq\\ [cos(a2;o — vr/4)Rej7'(m) — sin(a2;o — vr/4)Imj7'(ia)] , 

V t^zq 

where it follows from Eq. ( lA.llI) of Appendix A that 

Jiia) = -- e-^^"r(-za) . (39) 
a 



Comparison with Eq. f l37|) and the similar expression for Ji shows that to 
leading order, this result can be expressed in terms of Bessel functions 



A{s,zo) ~ zo [Jo{azo)ReJ'{ia) — Ji{azo)lmJ'{ia)] 



- + Im [e*^"r(m)] ) Ji{azo)+{Re [e'^'^T (ia)]) Jo{azo] 



a 



.(40) 



iii) Uniform approximation 

Expanding the Gamma function for small a, one obtains Eq. ( 136|) . showing 
that Eq. fHOl) gives a uniform approximation including both regions of the 
parameter a. Then from Eq. (!27|) it follows that for < a <^ l^ol, the 
asymptotic representation of the scattering amplitude is 

M{s,A) ~ { [[Im(e''^"r(2a))] J,{azo) + [Re{e'^''T (la))] Mazo)} . 

777/ -j^ 

(41) 
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5 Numerical results 



In this section we present some numerical results on the physical quantities 
of interest for different energy values. First let us come back to the deter- 
mination of the complex parameter zo{s) which plays an essential role in the 
solution of the high energy asymptotic behavior of the scattering amplitude. 
As we said earlier, zq is obtained by solving numerically Eq. f l30|) and the 
results are shown on Fig. 1, for Iiezo{s) and lmzo{s), versus Ins. As ex- 
pected, Rezo{s) rises rapidly with Ins, whereas Imzo(s) remains small and 
almost energy independent. It is worth noting that the asymptotic regime 
requires the validity of Eq. f l2^ . for example for Ins = lOvr, corresponding 
to the center-of energy y/s of about 6000TeV. 

The asymptotic representation of the forward scattering amplitude is ob- 
tained from Eq. fHT]) by taking the limit as a — )• 

v2 



Al(s,0)^3. (42) 



ISZq 

From this formula we calculate the ratio of the real to the imaginary parts 
of the forward amplitude defined as 



Re.M(s,0) 

P^'^ = fa^A^M) • 

In Fig. 2 (top) we display this result compared to the exact BSW result. We 
see that p(s) decreases for increasing energy, in agreement with the expecta- 
tion that p(s) — )■ 0, when s goes to infinity. 

The total cross section is obtained from the optical theorem as follows, 

a,ot = —lmM{s,0), (44) 
s 

and we recall that A^(s,0) is dimensionless. It is plotted in Fig. 3 (top) 
compared to the exact BSW result. 

In Figs. 2,3 (top), a gap can be noticed between the asymptotic representa- 
tion and the BSW model. The gaps extend to the end of the plotted energy 
range of s^^^ = 10^ TeV, where Ins ~ 37. To show that the asymptotic 
representation actually approaches BSW at sufficiently large values of In s, 
we carried out the first three terms of an asymptotic expansion to obtain 

is f 2 



Mis, 0) ~ ^ (^^o' + Y ) + 0(1/ zo) . (45) 
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With this expression in place of (42), the gaps largely close, as seen in Figs. 2,3 
(bottom). This shows indeed that the gaps visible in Figs. 2,3 (top) are due 
to the neglect of non-leading terms in the asymptotic representation. In Fig. 4 
we display the ratio of the leading order of the asymptotic representation to 
the exact BSW result, which goes to 1 for very, very large s, as expected. 
Now let us move from the forward direction to look at the behavior of the real 
and imaginary parts of the scattering amplitude as functions of t. For ^/s = 
14TeV, Fig. 5 displays the exact BSW amplitude along with its asymptotic 
representation. In both cases the imaginary part dominates the real part and 
its zeros will produce oscillations in the differential cross section, as shown 
in Fig. 6. The differential cross section is given by 

(46) 

where for the asymptotic representation one uses Eq. (HTI) . For both the 
BSW amplitude and its asymptotic representation, the real part has a lo- 
cal maximum near each zero of the imaginary part, and vice versa. When 
the maximum of the real part near a zero of the imaginary part is rela- 
tively low, as in the case near \t\ = O.SGeV^, one gets a sharp dip, but if 
not, like near \t\ = 2GeV^, one gets instead a smooth oscillation. Clearly 
the asymptotic result is larger than the exact BSW result, except near the 
diffraction peak, where they are hardly distinguishable. At a much higher 
energy ^/s = 6000TeV, the number of zeros increases, as shown in Fig. 7 
and the low maximum of the real part near the zero of the imaginary part 
around |t| = 0.2GeV^, generates a very sharp dip in the cross section, as seen 
in Fig. 8, followed by another dip and some smooth oscillations. 



6 Concluding remarks 

We have obtained the asymptotic representation of the BSW model in terms 
of a simple formula. The existence of several zeros for both the real and 
the imaginary parts of the scattering amplitude, generates oscillations in 
the differential cross section. The exact BSW result tends to coincide with 
this asymptotic representation, as the energy increases. This is even more 
striking near the forward direction, in particular for the diffraction peak of 
the differential cross section. In connection with this, let us mention now the 
following interesting feature of the asymptotic representation. A unitarity 
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upper bound for the imaginary part of the scattering amphtude for very 
high energy and small momentum transfers was derived a long time ago by 
Singh and Roy [7]. It reads Q 



lmMis,t) ^ 2Ji(v/f) 
ImA^(s,0) 



< 



if r < 3.46 



(47) 



with r = \t\atot{s) / 4:TT . For s very large, from Eq. ( 14T]) . we see that this vari- 



able is simply rj = aHezQ. The ratio 



lmM(s,t) 
lraM{s,0) 



has been plotted in Fig. 9, versus 



rj, for ^/s = 14TeV. We compare the exact BSW result with the asymptotic 
representation and also with the Singh- Roy upper bound limit s — )■ oo. We 
observe that the validity of the bound is, indeed, limited to i] < 3.46. Fig. 10 
displays the situation at a/s = 6000TeV and in this case the upper bound, 
whose validity is still limited to the diffraction peak, becomes much closer to 
the other two curves. 
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A Appendix 

For the small a region, we define 



dx 



-CXD 

■OO 



ln/3- 



d 

dv 



n poo 



dt fe 



u=0 



7 \ n 

ln/3--j T{l + u) 



u=0 



(A.l) 
(A.2) 
(A.3) 
(A.4) 
(A.5) 



where t = Pe ^ and F is the Gamma function. The results for n = 0, 1, 2, 3 
are given in Table 2. Here 7 ~ 0.5772 is the Euler's constant and ( is the 



"'^One should remember that the Singh- Roy amphtude is twice the BSW amphtude 
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n 





1 
2 
3 


1 

ln/3 + 7 

(ln/3 + 7)2 + 7r76 

(ln/3 + 7)^ + 4(ln/3 + 7) + 2C(3) 



Table 2: Values of X„(/3) 



Riemann Zeta function. The choice of Zq made in Eq. ( 130|) draws on the fact 
that 

X,(e-^) = -7 + 7 = 0, (A.6) 

hence Eq. ( 129|) follows as 







axx—e 
ax 



ax 



1-e 



dx 



1-e' 







dx e 



(A.7) 
(A.8) 



For the large a region one needs J{ia) and J{—iQ) where we define 



J{ia) 





1 

ia 
1 

ia 



1-e' 



_p-(^+7) 



c/xe*'' 

-1 + 
[-1 + e-*^"r(l - ia 



_p-(a;+7) 



dx 



(A.9) 
(A.IO) 
(A.ll) 



From Eq. f lA.lip and the properties of the Gamma function [6], it follows 
that J{—ia) is the complex conjugate of J{ia). 
For small values of a one finds 



Jiia) ^ C 3 . 



(A.12) 
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vs. In s 




Figure 1: The real and imaginary parts of Zq{s) as a function of the energy, 
obtained by solving numerically Eq. f l30|) . 
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p = Re[M(s,0)]/lm[M(s,0)] vs. s 



0.06^ 
10° 



BSW 

- - - asymp. rep 
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10* 



p = Re[M(s,C)l/lnn[M(s,0)] vs. s"' 



10 10 10 

s"^ [TeV] 



-BSW 

- asymp. rep. 



Figure 2: The ratio of the real to the imaginary parts of the forward ampU- 
tude versus y/s. Top using Eq. (H2l) and bottom using Eq. ( H5ll . 
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Ratio of asymp. rep. to BSW total scatt. cross sec vs. ^ 



0.! 



0.1 



0.' 




Figure 4: The ratio of the leading order of the asymptotic representation to 
the exact BSW result, versus the energy. 
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Figure 5: The absolute value of the real and imaginary parts of the elastic 
scattering amplitude, as a function of |t| for ^/s = 14TeV, for the exact BSW 
result and the asymptotic representation. 



20 




Figure 6: The elastic differential cross section calculated using Eq. ( I46i) . 
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p-p elastic scattering ampiitude for s = 6000 TeV 




Figure 7: The absolute value of the real and imaginary parts of the elastic 
scattering amplitude, as a function of |t| for a/s = 6000TeV, for the exact 
BSW result and the asymptotic representation. 



22 




Figure 8: The elastic differential cross section calculated using Eq. ( I46i) . 
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Normalized lm(Scatt. Amp) p-p elastic scattering amplitude for s = 14 TeV 



- imag(BSW) 

- imag(asymp. rep.) 
lim as s ^ «= 




n = Itl^'^Re z^m^ 



10 



12 



14 



Figure 9: The normalized imaginary part of the elastic scattering amplitude 
for = 14TeV. 
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Normalized lm(Scatt. Amp) p-p elastic scattering amplitude for s = 6000 TeV 
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Figure 10: The normalized imaginary part of the elastic scattering amplitude 
for = eOOOTeV. 
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